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TECHNICAL  NOTE  2982 

SUPERSONIC  FLOW  PAST  OSCILLATING  AIRFOILS 
INCLUDING  NONLINEAR  THICKNESS  EFFECTS 


By  Milton  D.  Van  Dyke 


SUMMARY 


A  solution  to  second  order  in  thickness  is  derived  for  harmonically 
oscillating  two-dimensional  airfoils  in  supersonic  flow.  For  slow  oscil¬ 
lations  of  an  arbitrary  profile,  the  result  is  found  as  a  series  includ¬ 
ing  the  third  power  of  frequency.  For  arbitrary  frequencies,  the  method 
of  solution  for  any  specific  profile  is  indicated,  and  the  explicit 
solution  derived  for  a  single  wedge . 

Nonlinear  thickness  effects  are  found  generally  to  reduce  the 
torsional  damping,  and  so  enlarge  the  range  of  Mach  numbers  within  which 
torsional  instability  is  possible.  This  destabilizing  effect  varies  only 
slightly  with  frequency  in  the  range  involved  in  dynamic  stability  analy¬ 
sis,  but  may  reverse  to  a  stabilizing  effect  at  high  flutter  frequencies. 
Comparison  with  a  previous  solution  exact  in  thickness  suggests  that 
nonlinear  effects  of  higher  than  second  order  are  practically  negligible. 

The  analysis  utilizes  a  smoothing  technique  which  replaces  the 
actual  problem  by  one  involving  no  kinked  streamlines.  This  strategem 
eliminates  all  consideration  of  shock  waves  from  the  analysis,  yet 
yields  the  correct  solution  for  problems  which  actually  contain  shock 
waves. 

I 


I NTRODUC  TI  ON 


As  linearized  supersonic -flow  theory  is  increasingly  applied  to 
problems  of  unsteady  motion  of  lifting  wings,  the  results  are  sometimes 
advanced  with  the  warning  that  they  may  be  significantly  affected  by 
nonlinear  effects  of  thickness.  Such  caution  is  justified  because  it  is 
known  that  even  for  steady  flow,  linearized  theory  is  often  inadequate 
for  predicting  the  pitching  moment  -  and  prediction  of  moments  is  one  of 
the  main  objectives  of  unsteady-flow  theory.  It  may  be  anticipated  that 
nonlinear  effects  will  become  increasingly  important  as  the  Mach  number 
falls  toward  unity,  particularly  for  slow  oscillations. 
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In  the  present  work  the  effects  of  thickness  are  determined  for  a 
harmonically  oscillating  two-dimensional  airfoil  by  calculating  the 
second-order  solution.  This  is  the  counterpart  for  unsteady  motion  of 
the  well-known  steady-flow  result  of  Busemann  (ref.  l) .  First,  for 
slow  oscillations  a  solution  is  found  for  an  airfoil  of  arbitrary  pro¬ 
file.  The  result  is  given  as  a  series  which  includes  terms  up  to  the 
third  power  of  the  frequency.  Second,  for  arbitrarily  high  frequencies 
it  is  shown  that  a  solution  can  be  found  for  any  specific  airfoil,  and 
the  solution  is  carried  out  explicitly  for  a  single  wedge.  Finally, 
comparison  is  made  with  a  previous  solution  for  the  wedge  which  is  exact 
with  respect  to  thickness  (refs.  2  and  3)>  in  order  to  assess  the  effects 
of  nonlinear  terms  of  higher  than  second  order. 

Extensive  use  is  made  of  a  smoothing  technique,  which  replaces  the 
actual  problem  by  one  having  no  kinked  streamlines.  This  strategem, 
which  has  been  used  previously  and  may  prove  useful  in  future  problems, 
eliminates  all  consideration  of  shock  waves  from  the  analysis.  It, 
nevertheless,  leads  to  the  correct  second-order  solution  for  the  actual 
problem,  which  does  involve  shock  waves. 


METHOD  OF  ANALYSIS 
Statement  of  Problem 


Consider  a  uniform  supersonic  stream  flowing  past  a  sharp-nosed 
airfoil  which  executes  prescribed  harmonic  oscillations.  We  shall  be 
concerned  with  calculating  the  instantaneous  pressure  at  the  surface  of 
the  airfoil  and,  hence,  the  unsteady  lift  and  pitching  moment.  It  is 
usually  the  custom  to  consider  the  oscillation  to  be  compounded  of  a 
rotation  (say  about  the  midchord)  and  a  vertical  translation  (plunging), 
which  are  not  generally  in  phase.  However,  we  shall  here  limit  atten¬ 
tion  to  the  special  case  of  the  airfoil  pivoted  about  an  arbitrary  point 
downstream  of  its  leading  edge.  This  is  no  real  restriction  because  the 
plunging  case  can  be  recovered  by  letting  the  distance  to  the  pivot  tend 
to  infinity  and  the  angular  amplitude  of  oscillation  tend  to  zero,  their 
product  remaining  finite. 

Although  the  iteration  procedure  to  be  employed  yields  a  formal 
result  for  any  Mach  number  greater  than  unity,  the  solution  probably 
breaks  down  when  the  flow  becomes  sonic  at  any  point.  Since  this  occurs 
at  a  Mach  number  somewhat  higher  than  that  for  bow  wave  detachment,  the 
upper  and  lower  surfaces  of  the  airfoil  operate  independently  in  the 
probable  range  of  validity  of  the  solution.  It  is  therefore  sufficient 
to  consider  only  the  half  field  of  flow  lying  above  the  airfoil,  and 
this  viewpoint  will  be  adopted  henceforth. 


NACA  TN  2982 


3 


Choose  the  origin  of  coordi¬ 
nates  at  the  mean  position  of  the 
leading  edge,  with  the  x  axis 
extending  in  the  direction  of  the 
free  stream  (sketch  (a)).  Then 
it  is  convenient  to  describe  the 
upper  surface  of  the  airfoil  in 
its  mean  (zero  angle  of  attack) 
position  by 


y 

4 


Sketch  (a) Coordinate  system  for 
oscillating  airfoil. 


y  =  Y(x)  =  e  g(x)  (l) 

Here  e  is  a  small  parameter  representative  of  the  airfoil  thickness, 
so  that  the  function  g  is  of  order  unity.  Now  let  the  airfoil  pivot 
about  a  point  lying  a  distance  b  downstream  from  the  leading  edge,  and 
perform  harmonic  oscillations  of  frequency  oj  and  amplitude  aQ,  so 
that  the  angle  of  attack,  which  is  the  angle  between  the  instantaneous 
position  of  the  airfoil  and  its  original  mean  position  (sketch  (a)),  is 
given  by 


a  =  aQ  cos  wt  =  (2) 

(Here,  as  in  all  that  follows,  it  is  implied  that  actual  physical 
quantities  are  given  by  the  real  parts  of  their  complex  representations.) 
Then,  at  any  instant  the  moving  upper  surface  of  the  airfoil  is  described 
implicitly,  to  second  order  in  thickness  and  angle  of  attack,  by 

y  =  e  g(x-ay)  -  (x-b)a  (3a) 

Now,  if  the  airfoil  is  smooth  (so  that  the  function  g  is  differenti¬ 
able),  the  equation  of  the  surface  can  be  written  explicitly  as 

y  =  e  g(x)  -  <xoeit0t  (x-b)  (3b) 

with  an  error  of  order  (e2a,ea2),  which  is  of  third  order  and,  conse¬ 
quently,  negligible  in  the  present  second-order  analysis. 


Perturbation  Equation 


The  entropy  changes  due  to  shock  waves  are  of  third  order  in  the 
airfoil  thickness  and  angle  of  attack.  Hence,  to  second  order  the  flow 
is  irrotational  and  isentropic.  Because  it  is  irrotational,  there  exists 
a  potential  function  ft  whose  gradient  yields  the  velocity  vector: 


k 
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T=  grad  ft,  v  =  fly  |  (4) 

Bernoulli Ts  equation  for  plane  unsteady  flow  can  be  written  (from 
eqs.  (l4.04)  and  (9*06)  of  ref.  4) 

ft+  +  i  (u2  +  V2)  +  —  =  +  1  u2  (5) 

2  7-1  7-1  2 

Here  a  is  the  speed  of  sound,  and  aQ  its  value  in  the  free  stream, 
where  the  flow  velocity  is  U.  Differentiating  this  expression  with 
respect  to  time  t,  and  using  the  fact  that  h[a2/(7-l)]  =  a2dp/p 
(ref.  4,  eqs.  (9*03)  and  (9.06)  gives 


fttt  +  uut  +  vvt  +  —  Pt  =  0  (6) 

This,  together  with  the  corresponding  results  obtained  by  differentiating 
with  respect  to  x  and  y,  can  be  used  to  eliminate  derivatives  of  the 
density  from  the  continuity  equation  (ref.  4,  eq.  (7*08.2)) 

Pt  +  (PU)X  +  (pV)y  =  0  (7) 

The  result  is  that  the  velocity  potential  satisfies  the  equation 

(a  -r2x2)Hxx  +  (^“^y2)^  "  ^^x^y^xy  "  ^^x^xt  ~  2Hyftyt  -  fttt  “  0  (8a) 

where,  from  equations  (4)  and  (5), 

a2  =  ao2  +  (U2  -  ftx2  -  fty2  “  2ftt)  (8b) 


Now  introduce  a  perturbation  potential  <J>,  normalized  through  division 
by  the  free-stream  velocity  U,  by  setting 


n  (x,y,t)  =  U  [x  +  $  (x,y,t)] 


(9a) 


so  that  the  velocity  components  are  given  by 


-  =  1  +  $ 
U 

v  & 

u  " 


x 


(9b) 
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Then  substituting  into  equations  (8)  gives 


(1-M2)  4^  +  cp^y-2  =  M2  (7-1)  (®x  +  ^  ~ (®xx+®yy)  + 


(2®X+4X  )  <^>XX  +  $yy  +  2  (  1+®X)  ^y^Xy  + 


I  <®x®xt  +  Vyt 


For  purposes  of  a  second-order  solution  (and  to  higher  order  a  potential 
does  not  exist),  the  triple  products  on  the  right-hand  side  can  he  dis¬ 
regarded.  Thus,  the  perturbation  equation  becomes  finally 


®  -P2®  ■ -  2 


yy  K  xx  U  ~xt  riStt 


vhere  p2  =  vF  -  1. 


®tt  =  M2  ( 7-1)  ( *v+  — )  ( 


U  xx1  -yy- 


2G>x<I>xx  +  S'Jy^xy  +  ^  (®x®xt  +  ^y^yt)  (ll) 


Pressure  Relation 


Dividing  the  Bernoulli  equation  (eq.  (5))  by  a?/{y-±)  gives 

_a£=  q  +  2zlM2A-  Bg+vg.+2fltN 

aQ2  2  \  U2  y 


The  flow  is  isentropic  to  second  order,  so  that 


Po  V  ao 


JPcU2 


and  it  follows  that  the  pressure  coefficient  at  any  point  is  given  by 


C  =  p~p°  =  JL 
15  1  p  U2  7M2 

2  0 


1  +  Z±  #  (l-  r'1-l 
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Substituting  the  velocity  components  of  equation  (9b),  expanding  in 
series,  and  retaining  only  squares  and  products  of  perturbation  quanti¬ 
ties  gives  finally 


CP  -  -  s*x  -2  Y  *  -•/  *  F  +  ^  (  r)2  d1*) 

Here  the  second-order  solution  is  required  only  for  evaluating  the  first 
two  terms;  the  others  are  given  correct  to  second  order  by  linear  theory. 


Smoothing  of  Problem 


From  the  leading  edge  and  from  any  subsequent  corner  of  the  airfoil 
spring  shock  waves  or  Prandtl-Meyer  expansion  fans  which  oscillate  as 
the  airfoil  oscillates.  These  introduce  serious  complications  into  the 
second-order  analysis.  However,  the  complications  can  all  be  circum¬ 
vented  by  solving  a  "smoothed"  problem  in  place  of  the  actual  problem. 
The  solution  can  thereupon  be  applied  to  the  actual  problem,  for  which 
it  yields  the  correct  result  everywhere  except  near  the  shock  waves  and 
Prandtl-Meyer  fans. 


Sketch  (b).-  Steady  flow  past  a 
wedge . 


The  nature  of  the  difficulties 
can  be  understood  by  considering 
first  the  special  problem  of  steady 
flow  past  a  single  wedge  (sketch  (b)). 
The  presence  of  the  bow  shock  wave 
means  that  the  analysis  must  be 
undesirably  complicated  by  including 
the  Rankine-Hugoniot  relations  (in  a 
simplified  form) •  A  second  compli¬ 
cation  arises  in  the  differential 
equation  which,  for  steady  flow, 
becomes 


$ 


yy 


-  P2®**  =  M2  (7+I) 


Mxx  + 


(r-l) 


®x®yy 


2®y^xy 


(15) 


In  the  iteration  procedure  to  be  employed,  the  nonlinear  right-hand  side 
is  evaluated  in  terms  of  the  first-order  solution,  and  the  resulting 
nonhomogeneous  wave  equation  solved  for  the  second-order  potential. 
However,  for  the  wedge  the  right-hand  side  vanishes  (to  any  order),  which 
would  imply  incorrectly  that  the  second-order  solution  does  not  involve 
the  adiabatic  exponent  7.  More  precisely,  the  right-hand  side  vanishes 
everywhere  except  along  the  Mach  lines  springing  from  the  apex,  where  it 
has  the  singular  behavior  of  the  Dirac  delta  function,  and  only  by 
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taking  account  of  these  troublesome  singularities  could  the  correct 
solution  be  found. 

Both  these  complications  are  avoided  by  the  simple  device  of 
solving  the  problem  of  flow  past  a  smooth  cusp-nosed  airfoil  of  arbi¬ 
trary  shape  and  then  applying  the  final  solution  to  the  wedge.  It  may 
be  imagined  that  the  wedge  has  been 
smoothed  by  adding  a  cusped  exten¬ 
sion  to  its  nose,  as  indicated  in 
sketch  (c).  It  is  clear  that  this 
artifice  removes  the  troublesome 
singularities  from  the  right-hand 
side  of  the  differential  equation. 

Likewise,  it  eliminates  the  need 
for  the  shock-wave  relations  be¬ 
cause,  as  indicated  in  sketch  (c), 
with  sufficient  smoothing,  shock 
waves  will  form  only  at  such  great 
distance  that  their  effects  cannot 
reach  the  airfoil  surface.  Although  Sketch  (c).-  Steady  flow  past 
shock  waves  are  thus  apparently  a  smoothed  wedge, 

excluded,  the  correct  second-order 

result  for  the  wedge  is  nevertheless  recovered  from  the  solution  by 
imagining  the  extension  to  shrink  in  size  and  disappear.  The  reason  is 
that  to  second  order  a  shock  wave  is  equivalent  to  the  limit  of  a  rapid 
continuous  isentropic  compression.  This  limiting  procedure,  which  is 
equivalent  simply  to  applying  the  solution  for  an  arbitrary  smooth  shape 
to  one  which  is  not  smooth,  yields  the  proper  result  except  in  the  vicin¬ 
ity  of  the  shock  wave  (see  ref.  5).  For  an  airfoil  of  general  shape, 
similar  broad  smoothing  must  be  imagined  at  any  concave  corner,*  whereas 
at  convex  corners  (since  no  shock  waves  form)  the  slightest  rounding  is 
enough.  This  smoothing  technique  was  applied  in  reference  4  (p.  399)  to 
steady  first-order  flow  past  bodies  of  revolution,  and  in  reference  5  to 
steady  second-order  plane  flow. 

We  turn  now  to  the  question  of  generalizing  this  smoothing  scheme 
to  an  oscillating  airfoil.  Modification  is  necessary  only  at  the  lead¬ 
ing  edge.  Consider  first  the  special  case  of  rotation  about  the  leading 
edge.  Then  it  is  enough  to  conceive  of  an  extension  which  is  flexible, 
so  that  its  cusped  tip  can  be  maintained  fixed  and  directed  always  into 
the  free  stream  while  the  airfoil 
oscillates,  as  indicated  in 
sketch  (d) .  (The  exact  motion  of 
the  flexible  tip  is  immaterial, 
provided  the  surface  is  suffi¬ 
ciently  smooth  and  its  slope 
remains  small.)  After  the  solution 
has  been  found,  the  flexible 
extension  is  again  imagined  to 


Sketch  (d).-  Smoothing  for  airfoil 
oscillating  about  leading  edge. 
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shrink  away,  and  the  correct  result  is  recovered  for  the  actual  airfoil 
oscillating  about  its  nose. 


Sketch  (e).-  Smoothing  for  airfoil 
oscillating  about  arbitrary 
point . 


Finally,  consider  rotation 
about  an  arbitrary  point.  The  flex¬ 
ible  extension  must  now  oscillate  in 
such  a  way  that  its  tip  is  alwayB 
directed  into  the  relative  wind. 
Hence,  as  indicated  in  sketch  (e), 
the  tip  must  lie  parallel  with  the 
free  stream  at  the  top  (and  bottom) 
of  each  stroke  but  incline  in  the 
direction  of  motion  for  intermediate 
positions . 


We  are  accordingly  led  to  consider  the  motion  of  an  arbitrary  flex¬ 
ible  oscillating  surface  described  by 

y  =  e  g(x)  -  aoeiu5t  h(x)  (l6a) 


where  for  the  smoothed  problem  the  functions  g(x)  and  h(x)  have  con¬ 
tinuous  first  derivatives.  The  smoothed  problem  will  ultimately  be 
replaced  by  the  actual  problem.  According  to  comparison  with  equa¬ 
tion  (3b),  this  means  that  the  function  h(x)  will  eventually  be  identi¬ 
fied  with  (x  -  b).  The  requirement  that  the  leading  edge  of  the  smoothed 
shape  be  always  parallel  to  the  relative  wind  may  be  written  as 


g'(0)  =  0 

h’(0)  =  iiob/u 


(16b) 


(These  last  conditions,  as  well  as  the  requirement  that  g'(x)  be  con¬ 
tinuous,  must  be  relaxed  in  recovering  the  solution  of  the  actual 
problem.) 


Boundary  Conditions 


The  boundary  condition  at  the  surface  of  the  airfoil  is  that  the 
normal  component  of  velocity  is  zero.  For  anysurface  described  by 
s(x,y,t)  =  0  moving  through  a  velocity  field  this  condition  means 
that  the  substantial  derivative  of  S  (i.e.,  its  time  rate  of  change 
for  an  observer  moving  with  the  fluid)  vanishes  at  the  surface  (see 
ref.  6),  so  that 


St  +  q  •  grad  S  =  0 


at  S  =  0 


(IT) 
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With  velocity  components  given  in  terms  of  Q  hy  equation  (9b),  and 
for  the  smoothed  surface  described  by  equation  (l6a),  this  tangency 
condition  becomes 

<£>y  =  (l+<I>x)(egT  -  ocoe1  h»)  -  ^  ccoe  h  at  y  =  eg  -  h  (l8) 

where  g  s  g(x),  etc.  It  is  convenient  to  refer  this  condition  to  the 
axis  y  =  0  by  expanding  in  Taylor  series.  Keeping  only  terms  of  second’ 
order  gives 

®y  -  (l  +  ^xHeg1 


(Here  $x  and 
ized  theory.) 

The  upstream  boundary  condition  requires  that  in  the  actual  problem, 
the  Rankine-Hugoniot  relations  (or  at  least  a  simplified  second-order 
form  thereof)  be  satisfied  across  an  oscillating  bow  shock  wave  whose 
position  must  be  determined.  However,  shock  waves  have  been  eliminated 
from  the  smoothed  problem,  so  that  it  is  only  necessary  to  require  that 
the  perturbation  potential  $  vanish  along  the  oscillating  character¬ 
istic  line  (Mach  line)  springing  from  the  leading  edge.  This  insures 
that  all  disturbances  produced  by  the  airfoil  are  swept  downstream. 

An  equivalent  and  still  simpler  requirement  is  that  $  and  its  stream- 
wise  derivative  vanish  on,  say,  the  plane  x  =  0: 

$=<j>x  =  Q  at  x  =  0  (20) 


„  _iwt  ico  icot  ,  /  iwt  _  \ 

"  Ooe  h')  -  —  a0e  h  -  (eg  -  c^e  h)$yy 

at  y  =  0  (19) 

on  the  right-hand  side  can  be  evaluated  from  linear- 


Restriction  to  Small  Angles  of  Attack 


It  is  convenient  to  seek  a  solution  to  second  order  in  the  airfoil 
thickness,  but  to  only  first  order  in  angle  of  attack.  This  is  not  a 
serious  limitation  because,  just  as  in  steady  flow,  second-order  terms 
in  angle  of  attack,  though  affecting  local  pressures,  have  no  effect 
upon  lift  or  moment  because  they  are  equal  on  the  upper  and  lower 
surfaces. 

As  a  consequence  of  this  restriction  to  linear  terms  in  angle  of 
attack,  the  awkward  question  of  when  the  bow  shock  wave  detaches  will 
not  arise.  The  amplitude  of  oscillation  may  be  regarded  as  infinitesi¬ 
mal,  so  that  if  the  bow  shock  is  attached  in  the  mean  steady  flow,  it 
will  remain  so  during  the  oscillation. 
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Transformation  of  Perturbation  Equation 

It  is  convenient  to  separate  the  time -dependent  part  of  the  problem 
from  the  mean  steady  flow  at  zero  angle  of  attack  (for  which  the  second- 
order  solution  is  known) .  Furthermore,  for  harmonic  oscillations  the 
number  of  independent  variables  is  then  reduced  to  two  by  separating  an 
exponential  time  factor.  Finally,  the  linear  portion  of  the  time- 
dependent  equation  is  reduced  to  normal  form  by  a  transformation  of 
dependent  variable.  These  three  transformations  amount  to  setting 

®(x,y,t)  =  <f>(x,y)  +  a0e1^tjt-KX^  (x,y)  (21a) 

where 

K 

p2  u 

Here  $  corresponds  to  the  mean  steady  flow,  and  the  term  in  $  repre- 
sents  the  additional  flow  associated  with  the  oscillation  through  small 
angle  of  attack.1 

Introducing  this  transformation  into  the  perturbation  equation 
(eq.  (ll))  gives  for 

*yy  -  P2*xx  =  M2  (7-D  <t>xA4>  +  (*x2  +  *y2)x  1  (22) 

where  A  is  the  Laplacian  operation  c^/dx2  +  c^/dy2,  and  for  $ 

$yy  "  P^xx  -  (if)  $  = 

ivr2  (7-l)(<l>xAf  +  $XA4>)  +  2(<t>^x  +  <t>y^y)x  -  iK^  (7-1)(24>x^x  +  ^j/A*)  + 
2  (“r1  +  <t>xx^+  ^  ]  ■  k2  [2  +  (7_1)  *  (23) 

1If  second-order  terms  in  angle  of  attack  were  retained,  equation  (21a) 

would  include  also  a  term  aQ2  e2i(wt-Kx)  x(x,y),  which  would  appear 

only  in  the  second-order  solution.  The  differential  equation  for  X 
would  have  the  form 

Xyy  -  P2Xxx  "  (t)  X  =  M2  *x  A  *  +  -  .  . 

and  could  be  solved  just  as  is  the  equation  for 


(21b) 
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The  tangency  condition  of  equation  (19)  likewise  separates  into  the 
two  conditions 

<t»y  =  e  (i  +  <t>x) g '  -  6  Ifryy  g  at  y  =  0  (24) 

=  -e1KX(l  +  4>x)  h'  -  iK  h-  e1KX  h  +  e(^x  -  iK$  g»  - 

e  ^yy  S  +  elKx  <t>yy  h  at  y  =  0  (25) 

For  the  actual  problem  the  second  of  these  becomes,  identifying  h(x) 
with  (x  -  b), 

fy  =  -eiKX  (l  +  <t>x)  -  iK  ^  elKx  (x  -  b)  +  e  (^-iieijr)  g'  - 

€$yy  S  +  elKX  'I’yy  (x  -  b)  at  y  =  0  (26) 


For  pressures  at  the  surface  of  the  airfoil,  the  relation  of  equation  (l4) 
can  be  expressed  in  terms  of  values  at  y  =  0  by  Taylor  series  expansion, 
with  the  result  that  to  second  order  in  thickness  and  first  order  in  angle 
of  attack 


CPu  =  (-24>x  -  2e<t’Xy  g  +  (324>x2  -  *y2)  + 


2a0ei(wt-ltx) 


P2*x$x  -  »y^  +  eiKX4>xy  h 


(27) 


Solution  by  Iteration 


Although  the  equation  for  4>  is  nonlinear,  that  for  ^  is  linear, 
hut  with  nonconstant  coefficients  depending  upon  <p.  This  corresponds 
to  the  physical  concept  that  because  of  the  restriction  to  linear  terms 
in  angle  of  attack,  the  oscillatory  part  of  the  flow  is  an  acoustic 
field  with,  however,  the  speed  of  sound  varying  from  point  to  point  in 
accordance  with  the  mean  steady  flow. 

The  well-known  linearized  or  first-order  theory  results  from  disregard¬ 
ing  the  right-hand  sides  of  equations  (22)  and  (23).  Thus,  with  the 
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first-order  potentials  denoted  "by  the  lower  case  letters  cp  and  i|r,  the 
perturbation  equations  become 

V  82  ’’xx  =  0  <28> 

*yy  -  P2  -  (if)  *  -  0  (29) 

The  second-order  solution  is  obtained  by  iterating  upon  the  first- 
order  results.  Using  the  linear  equations  to  simplify  the  right-hand 
sides  gives  for  the  second-order  iteration  equations 


^  yy  -  P2  *xx  =  2M2  [02  (W  -  l)cpx2  +  <p/] 


(30) 


and 

*yy  “  P2^ 


XX 


$  =  2M2  [  f  (N  -  1)  cpxtx  +  qyl^x  - 

/  \2 

21k  [02(2N-l)cpxtx  +P2NCPxx  t+cp  t  ]  -  2N  (£7) 

y  y  \m  /cp= 


2 

■* 

(31) 


where 


N 


=  7  +  1  W 
2 


(32) 


(Here,  following  the  usual  subscript  notation  for  derivatives,  [ty  2]  x 
means  5ly2/Sx,  etc.)  The  second-order  solution  for  <t>,  which  leads  to 
Busemann’s  well-known  result  at  the  airfoil  surface  (ref.  l),  was  given 
in  reference  5.  It  is  therefore  necessary  to  consider  only  the  second- 
order  problem  for  tJ.  Details  of  the  iteration  procedure  and  discussion 
of  its  limitations  are  given  for  the  steady  flow  in  reference  5  and 
apply  also  to  the  present  problem. 


Partial  Particular  Integral 


The  solution  of  the  differential  equation  for  steady  flow  in 
reference  5  was  simplified  by  discovery  of  a  particular  integral  of  the 
iteration  equation  in  terms  of  the  first-order  solution.  It  was  also 
shown  there  that  for  steady  three-dimensional  flow,  a  particular  integral 
can  be  found  to  account  for  all  terms  in  the  iteration  equation  except 
those  involving  the  adiabatic  exponent  y  in  the  form  of  N.  Likewise, 
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here,  a  partial  particular  integral  which  accounts  for  all  terms  on  the 
right-hand  side  of  equation  (31)  except  those  involving  N  is  given  by 

=  M2  (<Pt)x  "  iKCpt  (33) 

The  complete  solution  is  this  partial  particular  integral  plus  a 
solution  of  the  reduced  equation  whose  right-hand  side  contains  only 
the  terms  still  unaccounted  for: 

*yy  -  e%x  -  (f-J f  =  2p£?N  +  cpxtx]  -  ^  cpxt| 

(34) 


First-Order  Solution 


The  first-order  solution  for  cp  is  known  from  Ackeret's  theorv 

to  be 

?  =  “  J  g(x-Py)  (35) 


It  is  to  be  understood  here  and  in  all  similar  expressions  to  follow 
that  this  is  the  potential  only  for  x  >  Py,  and  that  cp  vanishes 
identically  ahead  of  the  bow  Mach  wave  (where  x  <L  Py) . 


The  first-order  equation  for  (eq.  (29))  is  most  readily  solved 
by  applying  the  Laplace  transformation  with  respect  to  x.  We  denote 
the  Laplace  transform  of  a  function  either  by  a  bar,  or  by  the  symbol  L 
whichever  is  more  convenient  (and  the  inverse  transform  by  L"1),  so  tha 


Applying  this  transformation  to  equation  (29),  using  the  fact  that 
and  i)fx  vanish  at  x  =  0,  gives 


(36) 

¥ 


ct  ^ 
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The  solution  of  this  equation  which  represents  waves  moving  downstream 
is 


t  =  c(s)  e-Py«/ s2  +  k2/m2 


The  constant  C(s)  is  determined  by  the  first-order  form  of  the  tangency 
condition  (eq.  (25)),  which  transforms  to 


ty  =  -  w(s)  at  y  =  0 


(39a) 


where 


w(x)  =  eiKX  I  h'(x)  +  in  pr  h(x) 


(39b) 


is  the  downwash  velocity  at  y  -  0. 
Consequently, 


_  1  _  e-Pyvs2  + 


l|f  =  —  w(s) 
p 


2  +  K2/M2 


The  inverse  transformation  is  readily  carried  out  using  the  standard 
tables  (e.g.,  ref.  7)  together  with  the  convolution  theorem,  which  gives 
as  the  solution  of  the  smoothed  problem 


x-Py 


i  (x,y)  =  j  J  J°  ^  J  v(x  -i)  d£ 


The  solution  for  the  actual  problem  is  now  obtained  by  setting 
h(x)  =  x  -  b,  which  gives  finally2 

x-Py 

*(x,y)=i  f  J„  ( eiK<x-»  |"l  +  Jk  £  (*-!-! 


iK(x-0  ["1  +  Ik  £  (x-i-b)l  d£  (42) 
M2  J 


in  agreement  with  the  known  result  of  linearized  theory  (see,  e.g.,  ref .  8) 


"The  smoothing  was,  of  course,  unnecessary  in  this  first-order  problem. 
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SECOND-ORDER  SOUJTION  FOR  LOW  FREQUENCIES 


Because  the  first-order  solution  for  arbitrary  frequencies  is 
rather  complicated,  use  is  sometimes  found  for  an  expansion  in  powers 
of  frequency,  which  involves  only  elementary  functions.  The  corre¬ 
sponding  second-order  solution  will  now  he  carried  out  in  detail, 
including  linear  terms  in  frequency.  This  result  will  serve,  for 
example,  to  evaluate  the  effects  of  thickness  upon  one -degree -of - 
freedom  torsional  instability,  which  is  primarily  a  low-frequency 
effect.  Thereafter,  the  result  of  extending  the  solution  to  include 
third  powers  of  frequency  will  simply  be  stated. 


Potential  Including  Linear  Terms  in  Frequency 


Expanding  the  first-order  solution  of  equation  (4l)  in  powers  of 
the  frequency  parameter  k  and  retaining  only  linear  terms  gives 

z 

*  =  J  f  w(i)  +  .  .  .  o(k2)  =  i|r(z)  +  .  .  .  (43) 


where  z  x  -  By.  To  this  order  the  partial  particular  integral  of 
equation  (33)  is  a  solution  of  the  homogeneous  equation  (eq.  (29))  and 
can  therefore  be  disregarded.  Substituting  the  first-order  solutions 
into  the  right-hand  side  of  equation  (34)  and  applying  the  Laplace 
transformation  gives,  to  order  k, 

^yy  "  p2s2f  =  -2eNe-Psy  (tfViiOL  jg'(x)\|r'(x)j  -Iks  L|g’(x)t(x)j  (44) 

It  is  readily  found  that  a  particular  integral  of  this  equation  is 
given  by  an  appropriate  multiple  of  ye-^.  Then,  adding  a  comple- 
mentary  function  representing  dovngoing  waves  gives 

where  the  constant  C  is>  to  be  evaluated  from  the  tangency  condition. 
Inverting  the  Laplace  transformation  shows  that  for  the  actual  problem, 
in  which  h(x)  =  x  -  b,  the  solution  has  the  form 


$  =  +  y  +  iK(2p2zg'  -  02bg'  -  g) 


(46a) 
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where  g  =  g(z).  The  arbitrary  function  f(z)  is  determined  from  the 
tangency  condition  of  equation  (26)  to  be 


f(z)  =  ^  +  ^ 


z  .  Ik  _  ( P 


p  P  \  2M2 


z  -  ^  b  +  e 


g  t  (x-b)  g' 

p 


±K€ 


2-(N-l)  2-^ 


,  nS  2^-1 

2(N_1)  ■  ww~ 


zg  +  z( z-b)g’ 


i  -  (N-b, 


(46b) 


Surface  Pressure  Coefficient 


The  pressure  coefficient  at  the  upper  surface  of  the  actual  airfoil 
is  found  from  equation  (27)  to  be 


Cpu  -  CPo  +  j-  1  +  *) 


M2N-2 


eg'  + 


IK 


(2-M2)(MsN-1) 

+  - S5 - xg’  +  (M2N-2)bg» 

P 


.  .  .  0( k2, e2a,a2) 


(^7) 


where  g  =  g(x).  Here  CpQ  is  the  value  for  the  mean  steady  flow  (at 
zero  angle  of  attack)  which  is  given  by  Busemann's  second -order  theory. 

A  more  useful  form  of  the  result  is  obtained  by  extracting  the  real  part 
and  expressing  the  result  in  terms  of  the  instantaneous  angle  of  attack 
a(t)  and  its  time  rate  of  change  a(t).  Furthermore,  the  parameter  e 
has  served  its  purpose  of  distinguishing  terms  of  different  orders  and 
can  be  eliminated  (according  to  eq.  (l)).  Thus,  on  the  upper  surface 
of  an  arbitrary  airfoil  which  is  described  at  zero  angle  of  attack  by 
y  =  Y(x),  is  pivoted  about  a  point  a  distance  b  downstream  of  its 
leading  edge,  and  performs  slow  angular  oscillations  described  by  a(t), 
the  pressure  coefficient  is,  to  second  order  in  thickness  and  first 
order  in  angle  of  attack. 
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2^(N-1)  y  +  (2-^)(MgIT-l)  xY,  +  t^N-2  by, 

P4  P4  p2 


Here  the  value  for  the  mean  steady  flow  is  (ref.  l) 


a 


~  +  •  ♦  .  0(a,e  a,ar 
(48a) 


2  ]ypN-2 
Cp  =  4  y>  +  Y’ 


o  P 


P‘ 


(48b) 


(A  preliminary  report  of  this  result  was  given  in  ref.  9*)  this 
form,  the  result  is  not  restricted  to  sinusoidal  motion  but  applies  to 
any  oscillation  which  is  sufficiently  smooth  and  slow  that  the  pressures 
depend  significantly  only  upon  the  instantaneous  angle  of  attack  and 
angular  velocity. 

The  pressure  on  the  lower  surface  of  the  airfoil  is  obtained  from 
these  equations  by  reversing  the  sign  of  a,  and  taking  Y(x)  to  be 
the  ordinate  of  the  lower  surface,  measured  positive  downward. 

The  result  for  plunging  motions  can  be  extracted  by  letting  a 
tend  to  zero  and  b  tend  to  infinity  in  such  a  way  that  their  product 
remains  finite,  say 

ba( t)  =  h(t)  (49) 


In  the  limit,  the  airfoil  simply  translates  vertically  according  to 
y  =  h(t).  The  pressure  coefficient  on  the  upper  surface  is3 


+  2£  +  2%2  Y* 
P  U  p2 


h 

U 


(50) 


Checks  on  the  Result 


The  solution  can  be  tested  in  several  special  cases  for  which  the 
result  can  be  derived  from  other  considerations. 


3  Note  that,  as  it  should  be,  this  is  just  the  result  of  using  Busemann's 
formula  for  steady  flow  (eq.  (48b)),  with  the  local  slope  increased  by 
the  instantaneous  apparent  downwash  angle  fi/U;  see  the  discussion  in 
the  following  section  of  conditions  at  the  nose. _ 
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Of  the  five  terms  in  equation  (48a),  the  first  is  known  from 
BusemannTs  steady  second-order  solution,  and  the  second  and  third  from 
linearized  unsteady  theory.  The  fourth  is  obtained  by  using  the  instan¬ 
taneous  airfoil  slope  (Yf-a)  instead  of  the  mean  steady  slope  Y1  in 
Busemann’s  formula  and  retaining  only  linear  terms  in  a.  Therefore, 
only  the  last  term,  which  is  the  essentially  new  result  of  the  present 
analysis,  requires  verification. 

Just  at  the  nose  of  an  oscillating  airfoil,  the  pressure  can  be 
determined  exactly  if  the  bow  shock  wave  is  attached.  The  transition 
through  the  moving  bow  shock  is  instantaneous,  and  so  depends  only  upon 
the  relative  velocity  at  that  instant  (see  ref.  4,  p.  297).  Hence  the 
pressure  just  at  the  nose  is  instantaneously  the  same  as  on  a  wedge  of 
the  same  vertex  angle  in  steady  flow  with  the  same  relative  velocity. 

In  the  present  problem,  the  relative 
velocity  is  compounded  of  the  hori¬ 
zontal  velocity  U  of  the  free 
stream,  plus  the  instantaneous 
vertical  velocity  of  the  leading 
edge,  which  is  given  by  db  (see 
sketch  (f)).  The  effect  of  the 
vertical  component  upon  the  equiva¬ 
lent  free -stream  velocity  and  Mach 
number  is  of  second  order  in  angle 
of  attack,  but  the  equivalent 
vertex  angle  of  the  airfoil  is 
increased  by  the  apparent  downwash 
angle  ab/u.  Replacing  Y*  by  Y*  +  ocb/U  in  Busemann*s  formula  (eq.  (Wb)) 
gives ,  to  first  order  in  angle  of  attack 


U 


Sketch  (f).-  Velocity  relative  to 
leading  edge  of  oscillating 
airfoil. 


!Pu  ~  CPo 


+  a  b  *  +  2  bY-  * 

p  u  p2  u 


(51) 


which  checks  the  part  proportional  to  b  of  the  last  term  in  equa¬ 
tion  (48a). 


Sketch  (g).-  Wedge  oscillating 
about  its  vertex. 


The  remainder  of  the  term  in 
question  can  be  checked  for  a  single - 
wedge  airfoil  oscillating  about  its 
vertex  (sketch  (g)).  It  can  be  shown 
using  the  results  of  reference  3 
that  in  this  case  disturbances 
reflected  from  the  shock  wave  are  of 
third  order  in  the  wedge  angle 
(although  for  other  pivot  positions 
they  are  of  second  order) .  There¬ 
fore,  a  solution  correct  to  second 
order  in  thickness  and  first  order 
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in  angle  of  attack  can  be  found  by  applying  linearized  theory  to  the 
mean  steady  flow  behind  the  shock  wave.4  For  slow  oscillations,  the 
first  three  terms  of  equation  (48a)  give 


|piUi2 


2  2(2-Mx2)  d 

-  a  +  -  Xi  — 

^  TT 


where  subscript  1  denotes  values  in  the  mean  steady  flow  behind  the 
shock  wave.  From  linearized  theory 


Mx  =  M  1-0  (N-l)  e 


Pi  =  p  [l  -  —  (N-l)  e] 

L  P  j 


Ui  =  U 


ri-ii 

pj 


Px  “  Po  1  +  7  —  e 

L  p 

r,  vf  1' 

Pi  =  Po  1  +  —  e 

Lei 


where  e  is  the  semivertex  angle.  Hence,  referring  the  surface  pressure 
coefficient  to  free-stream  quantities  (and  noting  that  Xi  =  x  to  second 
order)  gives 


n  _  n  2  „  ,  02 -M2  d 

LP„  =  -  -  a  +  2  — —  x  -  -  2 


M2N-2 


«  2Mb(N-1)  +  (2-MP) (M2 N-l)  A 

^  - ? -  €  X  ~ 

P  u 


which  checks  the  first  two  parts  of  the  last  term  in  equation  (48a) 
when  Y(x)  =  ex. 

Recently,  Lighthill  has  given  a  further  check  for  the  case  of 
Mach  numbers  so  high  that  l/l^  is  negligible  compared  with  unity 
(ref.  10). 

4 . . . —  "  - - - - - — - - - - 

This  concept  was  suggested  to  the  author  by  W.  P.  Jones  of  the  National 
Physical  Laboratory,  England. 
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Lift  and  Moment  Coefficients  for  Symmetrical  Airfoils 


The  coefficients  of  lift  and  pitching  moment  (about  the  pivot)  are 
given  in  terms  of  the  pressure  coefficients  on  the  upper  and  lower  sur¬ 
faces  by 

c 


(55a) 


cm 


(55b) 


For  simplicity,  consider  only  airfoils  symmetric  about  the  chord  line. 
In  this  case,  the  pressure  difference  is  given  by 


X  +  b  )  fi  +  4 
p2  x  +  D  j  D  +  4 


Y'a  - 


2Mg(N-l)  y  +  (2-Mg)(MeN-l)  xy,  +  t^N-2  byi 

-  P4  34  32 


a 


J  U 


(56) 


If  the  airfoil  has  a  sharp  trailing  edge,  substituting  into  equations 
(55)  and  integrating  by  parts  gives 


ie%li  54[T4l,2!f!!ii 

p3  c  Pc3 

o 


c 


Ydx 


o 


ccc 

U 


(58) 


Inside  each  square  bracket  the  first  term  is  the  result  of  linearized 
theory,  and  the  remainder  represents  the  second-order  effect  of  thick¬ 
ness.  Thickness  effects  are  seen  to  appear  in  the  form  of  the  area  of 
the  airfoil  profile  and  its  first  moment  about  the  vertical  line  through 
the  pivot. 
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If  the  airfoil  has  a  blunt  trailing  edge  of  semi thickness  Y(c), 
the  following  additional  second-order  terms  must  be  added  to  the  above 
expressions: 


to  c^: 


k  M2N-2  Y(c)  4 

3  3  c  a  p 


MfN-2  b  ,  (2-M2)(M2N-l) 

3  c 


Y(c)  ca 
c  U 


(59) 


to  cm: 

_  it  M2N-2 

3  3 


itWa.it 
V  e  3 


M2(N-1)  #N-2 


3s 


3 


^  1  A  Y(e)  cd 

c  J  j  \  cj  c  U 

(60) 


Examples:  Biconvex  and  Double-Wedge  Airfoils 


To  second  order,  a  biconvex  airfoil  of  thickness  ratio  t  is 
given  by  the  parabolas 

y  =  ±Y(x),  Y(x)  =  2t  |  (c-x)  (6l) 

The  expressions  for  lift  and  moment  become 


(63) 


It  happens  that  for  a  double -wedge  airfoil  (with  maximum  thickness  at 
midchord),  both  the  area  and  first  moment  are  just  three  fourths  of 
those  for  the  biconvex  airfoil.  Consequently,  the  above  results  apply 
to  double-wedge  airfoils,  if  t/3  is  replaced  throughout  by  r/k. 

In  the  expression  for  pitching  moment,  the  term  proportional  to  a 
represents  an  aerodynamic  stiffness  or  restoring  moment  in  phase  with 
the  angular  displacement,  while  the  term  proportional  to  d  corresponds 
to  an  aerodynamic  damping  moment  in  phase  with  the  angular  velocity. 

The  effects  of  thickness  upon  aerodynamic  restoring  and  damping  moments 
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are  shown  in  sketches  (h)  and  (i)  for  a  6-percent-thick 'double-wedge  or 
a  4-l/2-percent-thick  biconvex  airfoil  with  two  different  pivot  posi¬ 
tions.  The  sketches  have  been  labelled  with  both  the  usual  American 
notation  (e.g.,  ref.  ll)  and  the  British  notation  (e.g.,  ref.  12), 


LOO  L25  1.50  L75  2.00 


M 


<=>  __ ^ 


- Linearized  theory 

' —  Second- order  theory 


Jf 'lTeVt 

A  B  /  L  / 

Zi  N  G 

b/c  --  .5 

_ 1 _ 

. . . 

LOO  L25  L50  1.75  2.00 

M 


Sketch  (h).-  Effect  of  thickness  upon  restoring  moment  coefficient  for 
4-l/2-percent-thick  biconvex  or  6-percent-thick  double -wedge 
airfoil. 


- Linearized  theory 

-  Second-order  theory 


Sketch  (i).-  Effect  of  thickness  upon  damping  moment  coefficient  for 
4-l/2 -percent-thick  biconvex  or  6-percent  thick  double-wedge  air¬ 
foil. 


defined  for  harmonic  oscillations  of  arbitrary  frequency  by 


r 

-  |  k2  (Jfe  +  im4) 

(ma  + 


(American) 

\ 

(British) 

- 


(64) 
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Here  k  and  A  are  the  reduced  frequency  in  the  American  and  British 
notations,  respectively,  related  by 


A  =  2k  =  — 
U 


(65) 


For  slow  oscillations,  the  coefficients  are  given  by 

%  =  -  k2M3=  1 

2  da 


(and,  indeed,  this  was 
given  by  equation  (64) 
harmonic  oscillations) 


_U_  dcm 
2c  ^<5, 


(66) 


originally  the  definition  of  1%  and  m^,  that 
being  a  later  extension  to  the  case  of  rapid 


It  should  be  noted  that  according  to  second-order  theory,  the  non¬ 
linear  effects  of  thickness  are  themselves  linear  in  thickness.  This 
means,  for  example,  that  doubling  the  airfoil  thickness  ratio  would 
double  the  distance  between  the  linearized  and  second-order  curves  of 
sketches  (h)  and  (i). 


Neutral  Damping  Boundary 


Linearized  theory  indicates  the  possibility  of  instability  of 
pitching  oscillations  for  low  frequencies.  For  a  range  of  Mach  numbers 

below  ,/5/2  =  1.58  and  pivots  ahead  of  two  thirds  of  the  chord,  the 
aerodynamic  damping  moment  becomes  negative,  and  so  tends  to  destabilize. 
(Whether  or  not  the  motion  is  actually  unstable  depends,  of  course,  upon 
the  other  dynamic  parameters  in  the  problem.)  This  zone  of  possible 
instability  shrinks  and  eventually  disappears  as  the  frequency  of  oscil¬ 
lation  increases.  The  present  low-frequency  solution  is  therefore 
adequate  for  determining  how  the  region  of  instability  is  modified  by 
nonlinear  thickness  effects. 

Sketch  (j)  shows  the  boundary  of  neutral  aerodynamic  damping  for 
slow  oscillations  of  a  4-1/2 -percent-thick  biconvex  or  6-percent-thick 
double-wedge  airfoil.  The  aerodynamic  damping  is  destabilizing  for 
Mach  numbers  and  pivot  positions  lying  inside  the  loops.  Within  the 
region  where  linearized  theory  predicts  a  destabilizing  moment,  thick¬ 
ness  is  seen  to  exert  a  further  destabilizing  effect  except  for  pivots 
near  midchord.  The  second-order  solution  becomes  unreliable  when  the 
bow  shock  wave  detaches,  at  about  M  =  1.2. 
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Sketch  (j).-  Neutral  damping  boundary  for  4-l/2-percent-thick  biconvex 
or  6-percent- thick  double -wedge  airfoil. 


Comparison  With  Previous  Investigations 


Two  previous  investigators  have  sought  a  second-order  solution  for 
slowly  oscillating  airfoils  in  supersonic  flow.  Their  results  agree 
neither  with  each  other  nor  with  the  present  solution. 

In  1947,  W.  P.  Jones  obtained  an  estimate  of  the  thickness  effect 
by  assuming  that  the  ratio  of  second-order  to  linearized  pressure 
disturbances  is  the  same  for  slow  oscillations  as  that  given  by 
BusemannTs  formula  for  steady  flow  (ref.  13) •  That  this  assumption  is 
not  altogether  correct  is  indicated  by  the  fact  that  the  results  do  not 
check  those  obtained  for  a  wedge  oscillating  about  its  vertex  by  apply¬ 
ing  linearized  theory  to  the  mean  steady  flow  behind  the  shock  wave. 
However,  the  assumption  is  correct  at  the  leading  edge,  and  also  (as 
noted  by  Lighthill  in  ref.  10)  in  the  limit  of  high  Mach  number.  It  is 
seen  in  sketch  (j)  that  this  estimate  fails  to  give  a  useful  prediction 
of  the  actual  effects  of  thickness,  except  for  pivots  near  midchord. 

In  1951,  Alexander  Wylly  attacked  the  problem  by  methods  similar 
to  those  used  here  (ref.  l4) .  Unfortunately,  it  appears  that  the 
smoothing  was  not  carried  out  with  sufficient  care;  as  a  consequence. 
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the  solution  satisfies  none  of  the  three  checks  discussed  previously. 
In  contrast  to  the  present  results,  the  effect  of  thickness  upon  aero¬ 
dynamic  damping  was  predicted  to  he  stabilizing  and  so  great  that  for 
airfoils  of  the  thicknesses  shown  in  sketch  (j),  the  zone  of  possible 
instability  would  have  disappeared  altogether. 


Comparison  With  Experiment 


The  only  experimental  results  available  are  those  quoted  in 
reference  15,  which  are  derived  from  tests  by  Bratt  and  Chinneck. 

A  7-l/2 -percent -thick  biconvex  airfoil,  pivoted  at  midchord,  was  tested 
at  Mach  numbers  of  1.28,  1.46,  and  I.52  over  a  range  of  frequency.  The 
aerodynamic  damping  was  measured  by  giving  the  airfoil  an  initial  dis¬ 
placement  and  measuring  the  logarithmic  decrement.  The  experimental 
data  are  shown  in  sketch  (k)  in  comparison  with  linearized  theory  and 


¥ 


Sketch  (k).-  Comparison  of  experimental  and 
theoretical  damping  moment  coefficients  for 
7-l/2-percent-thick  biconvex  airfoil  pivoted 
at  midchord. 

the  second-order  theories  of  W.  P.  Jones,  Wylly,  and  the  present 
analysis.  As  indicated  on  the  sketch,  the  lowest  test  Mach  number  lies 
in  the  range  of  mixed  subsonic  and  supersonic  flow  (and  a  detached  bow 
wave).  The  present  theory  predicts  only  a  fraction  of  the  discrepancy 
between  experiment  and  linearized  theory.  The  same  is  true  of  Jones’ 
estimate  which,  according  to  sketch  (j),  agrees  with  the  present,  analysis 
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on ±y  for  pivots  near  midchord.  WyllyVs  result,  which  has  been  shown  to 
be  erroneous,  happens  to  lie  much  closer  to  the  measured  results. 


Extension  to  Cube  of  Frequency 


The  dependence  of  thickness  effects  upon  frequency  of  oscillation 
can  be  estimated  by  extending  the  second-order  solution  to  include 
higher  powers  of  frequency.  This  has  been  carried  out  for  an  arbitrary 
airfoil  by  including  second  and  third  powers,  which  is  enough  to  show 
an  effect  of  frequency  upon  both  aerodynamic  stiffness  and  damping. 

The  computation,  though  cumbersome,  is  a  straightforward  extension  of 
the  previous  analysis,  so  that  only  the  final  result  will  be  given  here. 


The  expression  for  pressure  coefficient  on  the  upper  surface, 
corresponding  to  equation  (47)  >  is  found  to  be 
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(67) 


X  X 

where  g  5  g(x)  and  / g  =  /g(l)d|,  etc.  (The  first  four  terms  are 

o  o 

the  result  of  linearized  theory.)  This  result  meets  the  three  tests 
discussed  previously,  and  also  checks  the  solution  given  later  for  a 
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wedge  with  general  pivot  location  at  arbitrary  frequency.  The  result¬ 
ing  expressions  for  lift  and  moment  involve  the  airfoil  thickness  in 
the  form  of  the  area  of  the  profile  and  its  first  three  moments  about 
the  vertical  line  through  the  pivot. 


Example:  Biconvex  Airfoil 


These  rather  formidable  results  simplify  considerably  for  specific 
airfoils*  For  example,  for  a  biconvex  airfoil  of  thickness  ratio  t 
oscillating  about  midchord,  the  pitching-moment  coefficient  is  given  by 


4a0elW't 


3  ^(6^-1 


l60j37 


^n-2  .  iAt  tflNrl).  + 
3P2  6P4 


A2tj^  12(2-M2)+(nP-4)N  +  lA3Tt^  2 (48m4 -61  ^+l)-(37M4+2^Mg-l6)N 
2^0 P 6  l440p8  lCQs 


The  first  three  terms  are  the  result  of  linearized  theory* 


The  component  of  this  moment  which  is  out  of  phase  with  the  angle 
of  attack  gives  a  parabolic  approximation  for  the  variation  of  aero¬ 
dynamic  damping  with  frequency  as  shown  by  dashed  lines  in  sketch  (z) 


Reduced  frequency ,  A  -  2k  -  uc/u 

Sketch  (Z)*-  Effect  of  frequency  upon  damping  moment  coefficient 
for  5 -percent -thick  biconvex  airfoil  pivoted  at  midchord. 
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for  a  5-percent-thick  biconvex  airfoil.  The  accuracy  of  the  parabolic 
approximation  for  linearized  theory  is  indicated  by  comparing  it  with 
the  exact  result  (solid  line) •  In  this  example  the  linearized  and 
second-order  curves  run  almost  parallel,  which  means  that  the  nonlinear 
effects  of  thickness  vary  only  slightly  with  frequency. 

Recently,  Jones  and  Skan  have  treated  biconvex  airfoils  at  arbi¬ 
trary  frequency  by  a  numerical  procedure  (ref.  16) .  Their  result  is 
shown  in  sketch  (z).  It  fails  to  give  the  initially  parabolic  form  of 
the  curve  which  is  implied  by  the  fact  that  the  second-order  solution, 
like  the  first-order  result,  can  be  expanded  in  powers  of  the  square  of 
the  frequency.  Their  solution  involves  several  doubtful  assumptions,  in 
particular,  that  the  effect  of  the  bow  shock  wave  can  be  disregarded. 

It  has  already  been  remarked  here  that,  actually,  the  bow  shock  has  a 
second -order  effect  unless  the  pivot  lies  at  the  leading  edge . 


SECOND-ORDER  SOLUTION  FOR  ARBITRARY  FREQUENCIES 


For  some  purposes  the  previous  solution  for  slow  oscillations  may 
be  inadequate.  In  principle,  the  second-order  solution  can  be  extended 
to  include  still  higher  powers  of  frequency,  but  the  labor  required  is 
clearly  prohibitive.  Alternatively,  one  can  attack  directly  the  problem 
for  arbitrary  frequencies. 

The  second-order  solution  can,  in  fact,  be  carried  out  for  a  general 
airfoil  at  arbitrary  frequencies.  However,  the  result  is  formidable, 
involving  multiple  integrals  of  products  of  Bessel  functions,  and  the 
reduction  to  simpler  form  for  specific  profiles  appears  to  be  difficult. 

A  more  practical  approach  is  to  choose  a  specific  airfoil  shape  in 
advance.  Then  the  second-order  solution  involves  only  functions  of  the 
type  encountered  in  the  linearized  theory.  In  particular,  it  is  found 
that  (at  least  for  the  simplest  shape)  the  final  expressions  for  lift 
and  moment  involve  only  functions  which  have  been  already  studied  and 
tabulated,  so  that  numerical  results  are  readily  obtained. 


Modified  Smoothing  Procedure 


The  smoothing  discussed  previously  must  be  dropped  at  an  earlier 
stage  of  the  solution  when  a  specific  airfoil  shape  is  chosen.  It  is 
therefore  necessary  first  to  modify  the  differential  equation  and 
boundary  conditions  so  that  no  singular  terms  appear. 

Consider  first  the  differential  equation  (eq.  (31))*  Applying  the 
Laplace  transformation  of  equation  (36),  and  envisioning  the  smoothed 
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problem,  so  that  $  and  $x  vanish  at  x  =  0,  reduces  it  to 
fyy  -f  f  *  2(mPs-1k)  |^P2(N-1)L-|9jC'I'x|  +  L^Pytyj--iK  ^ 

(69) 


In  this  form,  all  troublesome  second  derivatives  have  disappeared  from 
the  right-hand  side,  so  that  the  smoothing  can  now  be  dropped. 

Consider  next  the  tangency  condition  of  equation  (25),  which 
contains  second  derivatives  of  both  9  and  These  cause  no  difficulty 
when,  as  in  the  solution  for  slow  oscillations,  the  tangency  condition 
is  imposed  as  it  stands.  However,  its  Laplace  transform  will  be  used 
here,  and  then  the  Dirac  delta  functions  associated  with  the  second 
derivatives  would  affect  the  integration  implied  in  the  inversion  of  the 
transformation.  These  troublesome  second  derivatives  can  be  eliminated 
by  first  expressing  them  in  terms  of  x  derivatives  only  through  the 
first-order  equations: 


fyy  =  P2  %cx 

^yy  =  P  ^xx  +  JJp  ^ 


(70) 


and  then  applying  the  Laplace  transformation  while  envisioning  the 
smoothed  problem  (so  that  <p,  i,  cp*,  and  \|rx  vanish  at  x  =  0),  which 
gives 


-  w(s)  - 


Again,  all  second  derivatives  have  disappeared,  so  that  the  smoothing 
can  be  dropped. 

Consider  finally  the  upstream  condition.  For  the  smoothed  problem, 
one  statement  of  the  condition  was  seen  to  be  that  the  solution  repre¬ 
sents  downgoing  waves.  This  means  that  the  complementary  function  for 
the  iteration  equation  should  have  the  same  form  as  the  linearized 
solution,  and  this  statement  of  the  upstream  condition  applies  as  well 
to  the  actual  problem. 
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Solution  for  Wedge 


The  simple  case  of  a  single  wedge  illustrates  the  method  of  solu¬ 
tion  which  can,  in  principle,  be  applied  to  any  profile  formed  of 
piecewise  analytic  arcs.  For  a  wedge  of  semi vertex  angle  e,  the 
first-order  solution  for  <p  is  given  by  equation  (35)  as 

<p  =  -  £  (x  -  py)  (72 

P 


and  the  first-order  solution  for  t  is  given  by  equation  (42).  Substi 
tuting  into  the  right-hand  side  of  equation  (69)  gives  for  the  trans¬ 
formed  iteration  equation 


fyy  -P2  (s2+  -  ^(M2 


s-iK)  [l. 

L  J  s2+k2/M2 


w(s)e 


|3y  ,</s2+k2  /M2 


A  particular  integral  of  this  equation  is  given  by  an  appropriate 

multiple  of  ye-^^8  +K  y  and  adding  a  complementary  function  which 
represents  downgoing  waves  gives  the  solution 

C(s)e-Py^32+K2/M2  + 


-  (l^s-iic) 
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i2+K2/M2 


— (a)ye-Py»/  s2+k2/M2 
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The  constant  C  is  evaluated  by  imposing  the  tangency  condition  of 
equation  (71),  with  the  result  that  the  Laplace  transform  of  the  second- 
order  solution  is  found  to  be 


V  e-(3y  v  32+k2/M2  f-/  r  7  o  Go 

$  =  . - . -  ■  - w(s)  l+ey(l^s-iK)  +  -  (M2s-1k) 

pys2^2/^  1  l  p 


(N-l^-iKN/M2 


■2+k2/M2 


s^/m2 


_  Pe +k2/M2  A  w(s)  _  Peb  _  f  — i- 
ds  P  s-in 


^n  this  simple  case,  the  partial  particular  integral  of  equation  (33) 
offers  little  actual  simplification. 
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The  inversion  ean  be  carried  out  using  the  standard  tables  (e.g., 
ref.  7)  together  with  the  convolution  theorem.  For  calculating  surface 
pressures,  it  suffices  to  obtain  the  solution  in  the  plane  v  =  0. 
which  is  found  to  be 


if  (x,°)  =  J  jQ  ^  £  1^)  w(x-£)d£  +  ex  w(x) 
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(76a) 


where 


w(x)  =  eiKX 


1  +  iK 


If  <x-b) 


(76b) 


With  the  pivot  at  the  nose  (b  =0),  this  agrees  with  the  result  of  apply¬ 
ing  linearized  theory  to  the  mean  steady  flow  behind  the  shock  wave. 

Also,  when  expanded  in  powers  of  frequency,  it  agrees  with  the  previous 
low-frequency  solution  up  to  terms  in  K3. 

The  surface  pressure  coefficient  can  now  be  calculated  from  equa¬ 
tion  (27)  and  the  lift  and  moment  coefficients  from  equations  (55). 


Example:  Wedge  Pivoted  at  Nose 


For  simplicity,  the  results  will  be  given  only  for  the  special 
case  of  rotation  about  the  nose.  Then  it  is  found  that  the  lift  and 
moment  coefficients  are  given  by 
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(77) 


where  V  =  0  for  the  lift  and  v  =  1  for  the  moment.  The  functions 
fn>  given  hy 


fn(M,A) 


J  xHe-iA^x/p2  Jq 


(78) 


arise  in  linearized  theory  for  n  ranging  from  0  to  3*  They  have  been 
studied  and  tabulated  by  von  Borbely  (ref.  17) >  Schwarz  (ref.  18),  and 
Garrick  and  Rubinow  (ref.  ll) .  They  can  all  be  expressed  in  terms  of 
f0  by  a  recurrence  relation  due  to  von  Borbely  (see,  e.g.,  eq.  (A. 87) 
of  ref.  12),  so  that  the  additional  f4  required  here  is  easily 
computed. 
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Sketch  (m)  shows  the  variation  of  aerodynamic  damping  moment. 


O  A  . 8  L2  L6  2.0  2.4  2.8 

Reduced  frequency ,  A -2k-  »c/u 


Sketch  (m).-  Damping  moment  coefficient  for  10-percent- 
thick  wedge  oscillating  about  its  vertex  at 
M  =  10/7. 


according  to  first-  and  second-order  theories,  for  a  10-percent-thick 
single  wedge  oscillating  about  its  nose  at  a  Mach  number  of  10/7«  Also 
shown  for  comparison  are  the  parabolic  approximations  of  the  low-frequency 
analysis*  It  is  seen  that  the  thickness  effect  is  reversed  at  high 
flutter  frequencies,  as  is  suggested  by  the  parabolic  approximation. 


DISCUSSION 
Higher-Order  Effects 


The  moderate  magnitude  of  second-order  effects  would  suggest  that 
the  influence  of  third-  and  higher-order  terms  is  of  no  practical 
importance,  except  perhaps  in  the  transonic  range  near  shock  detachment. 
This  supposition  can  be  confirmed  in  the  case  of  the  single-wedge  airfoil, 
for  which  a  solution  exact  in  thickness  (but  linearized  with  respect  to 
angle  of  attack)  has  been  derived  in  references  2  and  3*  Sketch  (n) 
compares  the  boundaries  of  neutral  aerodynamic  damping  for  a  slowly 
oscillating  wedge  of  5°  semivertex  angle  as  predicted  by  the  linearized. 
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Pivot  position ,  b/c 


Sketch  (n).-  Neutral  damping  boundary  for  wedge  of  5° 

semivertex  angle. 


second-order,  and  exact  theories.  The  second -order  solution  lies  close 
to  the  exact  result  down  to  the  Mach  number  for  shock  detachment  (which 
is  almost  the  same  as  the  Mach  number  at  which  the  flow  ceases  to  be 
purely  supersonic). 


Application  to  Finite-Span  Wings 


Extension  of  the  second-order  solution  to  wings  of  finite  aspect 
ratio  does  not  seem  possible  at  present.  No  second-order  solution  has 
yet  been  found  even  for  steady  flow  past  the  simplest  lifting  wing. 

Fortunately,  the  main  conclusion  to  be  drawn  from  the  present 
analysis  is  that  nonlinear  thickness  effects  are  quite  moderate  in 
magnitude.  Practical  supersonic  wings  will,  therefore,  probably  be 
so  thin  that  nonlinear  effects  are  negligible,  so  that  reliance  can  be 
placed  in  the  predictions  of  linearized  theory.  Only  if  the  wing  is 
unduly  thick,  or  if  the  Mach  number  is  close  to  unity,  or  if  unusual 
accuracy  is  required,  may  the  engineer  be  forced  to  estimate  the  effects 
of  thickness.  In  this  event,  he  might  assume  that  the  effects  of  thick¬ 
ness  are  in  some  sense  additive  to  those  of  aspect  ratio,  provided  the 
aspect  ratio  is  high  and  the  frequency  low.  For  example,  the  two- 
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dimensional  correction  might  he  applied  stripwise  to  the  spanwise  load¬ 
ing  predicted  by  linearized  theory.  Some  indication  of  the  extent  to 
which  such  an  assumption  would  be  justified  can  be  obtained  by  consider¬ 
ing  other  pairs  of  effects  whose  combined  influence  is  known.  Sketches 
(z)  and  (m)  show  that  the  effects  of  thickness  and  frequency  are  roughly 
additive  for  the  frequencies  of  usual  practical  interest  in  dynamic 
stability  analysis  (say  A<  0.2) .  Likewise,  sketch  (o)  shows  that  the 


Sketch  (o).-  Effect  of  aspect  ratio  upon  damping  moment 
coefficient  for  rectangular  wing  according  to 
linearized  theory. 


effects  of  aspect  ratio  and  frequency,  determined  from  Watkins'  linear¬ 
ized  solution  for  the  rectangular  wing  (ref.  19),  are  nearly  additive 
in  the  same  range  of  frequencies. 

Recently,  Martin  and  Gerber  have  calculated  the  second-order  effects 
of  thickness  upon  damping  in  roll  for  wings  of  infinite  span  (ref.  20). 
They  then  estimate  the  effect  for  a  finite  rectangular  wing  by  increasing 
the  result  of  linearized  theory  in  the  same  ratio  as  for  the  infinite 
wing.  Agreement  with  experiment  is  thus  considerably  improved,  which 
gives  further  assurance  that  superposition  of  thickness  effects  may  be 
justified  also  for  oscillating  wings. 
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Further  Analysis 


If  thickness  effects  of  the  magnitude  indicated  by  the  present 
analysis  are  judged  to  be  significant  at  flutter  frequencies,  extension 
of  the  high-frequency  solution  to  more  practical  profiles  would  be 
warranted.  The  solution  given  for  a  single  wedge  at  arbitrary  frequency 
should  be  extended  next  to  the  biconvex  or  double -wedge  airfoil. 

Although  considerable  computation  is  involved,  it  does  not  appear  that 
the  labor  would  be  prohibitive. 


Ames  Aeronautical  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Moffett  Field,  Calif.,  April  20,  1953 


NACA  TN  2982 


37 


APPENDIX 

SYMBOLS 


a 

A 

b 

c 

C(s) 

cl 


speed  of  sound 
aspect  ratio  of  wing 

downstream  distance  from  leading  edge  to  pivot 

airfoil  chord 

constant  of  integration 

section  lift  coefficient 


CPU'  CPZ 
g(x) 


h(t)' 

h(x) 


k 


section  moment  coefficient 

pressure  coefficient 

surface  pressure  coefficient  in  mean  steady  flow 

pressure  coefficients  on  upper  and  lower  surfaces  of 
airfoil,  respectively 

function  defining  upper  surface  of  airfoil  at  zero  angle  of 
attack 

elevation  of  plunging  airfoil 

function  defining  amplitude  of  oscillation  of  upper  surface 
of  smooth  airfoil 

wc 

reduced  frequency,  — 

2U 

aerodynamic  damping  coefficients 
(See  eq.  (66)) 


M  free- stream  Mach  number 

n  Zii!£ 

2  g 

p  static  pressure 

q  velocity  vector 

s  Laplace  transformation  variable 
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S(x,y,t)  function  defining  moving  surface 
t  time 

u,v  velocity  components  parallel  and  perpendicular  to  free  stream 

U  free-stream  velocity 

w(x)  downwash  velocity  at  y  =  0 

x,y  coordinates  parallel  and  perpendicular  to  free  stream 

Y(x)  ordinate  of  upper  surface  of  airfoil  at  zero  angle  of  attack 

z  x  -  Py 

a  angle  of  attack 

ccq  angular  amplitude  of  oscillation 

0  a/m2  "  1 

7  adiabatic  exponent  of  gas 

A  Laplacian  operator,  -  +  - 

dx2  dy2 

e  small  parameter  representative  of  airfoil  thickness 


7  reduced  frequency,  wc 

U 

V  constant  vhich  is  0  or  1  in  equation  (74) 

p  density 

t  airfoil  thickness  ratio 

cp  first-order  mean  steady  perturbation  potential 

second-order  mean  steady  perturbation  potential 
$  complete  perturbation  potential 

t  first-order  time -dependent  perturbation  potential 

^  second-order  time -dependent  perturbation  potential 
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fl 


<  >o 

(  )i 
(  *) 
C) 


L 


partial  particular  integral  of  time-dependent  iteration  equation 
angular  frequency  of  oscillation 
complete  velocity  potential 
free- stream  conditions 

conditions  in  mean  steady  flow  behind  bow  shock  wave  on  wedge 
differentiation  with  respect  to  time 

Laplace  transform 
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